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Abstract 

We extend the model for QCD string with quarks to consider the Mercedes Benz 
string configuration describing the three-quark baryon. Under the assumption of 
adiabatic separation of quark and string junction motion we formulate and solve 
the classical equation of motion for the junction. We dare to quantize the motion of 
the junction, and discuss the impact of these modes on the baryon spectra. 

The qqq baryon, together with the qq meson, is usually considered in the framework 
of the constituent model as a simplest system with zero triality. It is clear, nevertheless, 
that in the QCD -motivated approach the only way to form gauge-invariant states with 
three quarks in the fundamental colour representation is to introduce the so-called string 
junction as an object that couples three quarks to make a colour singlet: 

^(3:1X2X3, Xo) = 

where parallel transporters $, 

^%x,y) = {Pexptg r Blt^dz^Tp , (2) 

stand to assure non-local gauge invariance of the qqq system placed into non-trivial QCD 
background -B^. From this point of view it is obvious that the baryon should be treated 
as a simplest multiquark system containing, in comparison to the qq system, a new object 
with it's own dynamics. Standard way to deal with the junction dynamics is to assume 
that it always moves in such a way that the sum of distances between quarks and junc- 
tion is minimal [1,2]. This assumption is motivated by non-relativistic approach, where 
confinement manifests itself as linear potential between constituents. On the other hand, 
the idea that a junction might be responsible for certain kinds of baryonic excitations was 
already discussed in [2]. 
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In the present letter we relax the above-mentioned assumption on a junction following 
submissively the quarks. Actually our baryon is going to contain four constituents rather 
than three. Our starting point is the Vacuum Background Correlators method [3]. In 
the framework of this method the hadronic Green functions can be constructed starting 
from the QCD Lagrangian. The dynamics of interaction is defined by the averages of 
Wilson loop operators, and it may be shown that the area law asymptotic for Wilson loop 
average gives rise to the string-type interaction of the constituents. The model for the 
resulting straight-line qq string was considered in [4], and the hybrid mesonic excitations 
were analysed in [5]. 

To formulate the Hamiltonian approach to the problem one is to consider the Green 
function of a baryon, 

^(xixaxaxo, ?/i?/2y3yo) =< ^(xiXsXa, Xq) | ^iymys^Vo) >b , (3) 

where the brackets mean the averaging over background field configurations. To define 
the effective action the Feynman-Schwinger representation [6] of a Green function (3) is 
used (in the Euclidean space): 

/•oo roo roo r 

G{xiX2X3Xo,yiy2y3yo) = dsi ds2 / ds^ / DZ1DZ2DZ3X 

Jo Jo Jo J 

X exp(-/C)- <W >B (4) 

where 

i=l ^ ''^ 

boundary conditions are -2j(0) = yi , Zi{si) = Xi, and the contours Fj in the Wilson loop 
operator run over the trajectories of the quarks. Using the relation 

1 = ^^«.«.«3^"^"^"^*^Kr)<^>^:(r)$^3'(r) (5) 

(which is due to the unimodularity condition for the SU(3)), where F is an arbitrary open 
path (the same for all three $'s) connecting points xq and yo, we represent the Wilson 
loop in (4) as 

W = SpWi SpW2 SpW3 - Sp{WiW2) SpW3- 

- Sp{W3 Wi) SpW2 - Sp{W2 W3) SpWi+ (6) 

+Sp{Wi W2 W3) + Sp{W3 W2 Wi) 

where (Wi)'^ is the ordered exponent along the contour formed by the paths of i-th quark 
and F.Q 

The averaging over background was done using the cluster expansion method [3], 
generalized for the case of more than one Wilson loop as well as for a loop with self 

^In what follows T will be interpreted as the path of the junction. 
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intersection in [7] . The main result of this method is that under assumption of existence of 
finite gluonic correlation length Tg the generalized area law asymptotics may be obtained. 
For the Wilson loop (6) we write 

< W >B- j DzoeM-'^Si - aS2 - aS^) , (7) 

where a is the string tension and Si is the minimal area inside the contour Cj, and we 
integrate in (7) over the junction trajectories {2:0} in accordance with our intention to treat 
a junction as a degree of freedom. The area law (7) is held for the contours with average 
size much larger than Tg, and is violated only when the contours arc nearly embedded 
into the same plane. In what follows we neglect such special configurations. 

We would like to note here, that the standard approach to the junction motion cor- 
responds to taking the classical trajectory Zo{r) in the integral (7) or, equivalently, as- 
suming J2f^i Si = mm. The latter condition can be reduced in the potential model to the 
assumption that the sum of distances between quarks and junction is minimal. In the 
"string-type" language the representation (7) means that we are interested in the special 
kinds of string excitations, which are absent for simple qq string configurations, and reveal 
themselves in multiquark systems. 

Using the parametrization 

Zif,^{T,fi), ZOi^^{T,fo) 

and introducing new dynamical variables 

l^iir) = ^z,o{t) 0<t<T , 

ZSi 

we rewrite the Green function (4) as 

G = J DfiDf2Df3DfoDfiiDfi2Dfi3exp{-A) (8) 



with the effective action 

3/^,2 , -2 



(9) 



where the surfaces Si arc parametrized by the coordinates -iUj^, and tUj^ = , w'i^ = 

-Q^- This procedure of reduction the four-dimensional dynamics to the three-dimensional 
one (9) was first suggested in [4] to analyse the qq system. With the straight-line ansatz 
for the minimal surfaces 

Win. = ~ A) + ^0/.(t)A, 

we obtain the Lagrangian 

(10) 
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- 1 • • 

= — [a X ((1 - pi)fi + pifo)] , pi^Vi-ro 
Pi 

The Lagrangian (10) with the additional requirement J2i=iPi — "^^^ was discussed in 
detail in [8] with neglecting the angular velocities li under the square roots (string cor- 
rection), and the attempt to account the string correction was undertaken in [9]. 

Here we concentrate on the dynamics of junction, so we assume the quarks to be heavy 
enough (m^ ^ ^/a) to allow the adiabatic treatment of the problem. Following again the 
proceduce of [4] we introduce the auxiliary fields 1^1,1]^. 

333 

G = ^tlD'y^llDr]^ll DpiDrMDnDfoexp{-A[{u], {r;}]), (11) 
i=i i=i j=i 

+^ / dPi'^iVii^i - ^of + d(3iUirii{fi - fo, ^fo + (1 - A)^0 j 

The integration over the fields z/j and rji actually amounts to taking extremal values 
in Vi and rji. 

To make the algebra as simple as possible we consider the case of equal quark masses. 
Moreover, we take into account only the lowest symmetric hyperspherical harmonic [10] in 
the quark subsystem, an approximation which is proved to be very successful in few-body 
systems. It is quite clear that under all these assumptions one has 

A*! = A«2 = A«3 = A«, i^i = = i^3 = J^, vi^m^m^ V, 

and, after integrating over the fields p and rj we arrive to the relatively simple Hamiltonian 
in the center-of-mass frame (in the Minkowski space): 

H = 3m + ?^±^ + H, (12) 
+ + (13) 

Jo 

where the quark wave finction depends only on the hyperradius p, p^ = {fi — + 
{f-i — '''^^"^ )'^ ^ r = fo, n = ^, and p, q and Q are the momenta conjugated to the Jacobi 
coordinates fi — f2, fa — ^^li^ and f correspondingly. 

The adiabatic approximation means that one should find the eigenenergies En of the 
"fast" subsystem (13) as functions of p, and after that substitute En{p) as adiabatic 
potentials into the Hamiltonian (12) of "slow" subsystem.To perform this program we 
first found the classical equations of motion for the Hamiltonian (13) with the result ( we 
restrict ourselves with the spherically symmetric case of L = 0) : 



Jo z/ ' 



m, = / apup', m.uP' — 



4 



Eia) = ]-m,uj^p\a^ + + l t '^(^^^ (14) 
I I Jo 

{Qn) = mup^-^^-^, (15) 

where a is the constant of integration. 

It may be shown from equations (14), (15) that "freezing" of the junction degree of 
freedom actually means substitution the lowest classical energy £"(0) into the Hamiltonian 
(12). In this case integration over i> is trivial, and one arrives to the potential model 
Hamiltonian [1]0 

H = 3m + ^ ^ + VSap (16) 
3m 

On the contrary, the Bohr-Zommerfeld quantization procedure for the classical solution 
(14), (15) yields 

1 3 /"^ 

En = -m,tuV(a^ + 1) + - / df3iy, (17) 



where 



msUJp^J{an) = (2n+ ^)^, 



a 



la? — x'^ 



and may be expressed in terms of elliptic integrals K(-yj=5=) and E( ^^^^^ ). The integra- 
tion over u now cannot be performed analytically, apart from the case of asymptotically 
large p, p ^ where the spectrum (17) becomes the spectrum of harmonic oscillator, 
and the adiabatic potential takes the form 

En{p) = V3ap+{2n + \)--. (20) 

z p 

We can find, however, the value of p that corresponds to the minimum of the potential 
well for quarks: 

2 7r(2n+|) 

2m,cu2j(a(0))' ^ > 

where for any n a*^°^ ~ 2.2 is the root of the equation 

2aJ(a)-(a2 + 1)^^ = (22) 
oa 



^The coefficient at the confinement form in (16) differs siightly from one obtained by taking lowest 
hyperspherical harmonic of the interaction with a X]i=i Pi — ''min. It happens because we integrate over 
the field v after hyperspherical decomposition rather than before. The order of these procedures does 
not matter only if one does not truncate the basis. 
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Performing the integration over u for po defined from the equation (21), we arrive to 



ext 



Having in mind that our quarks are heavy( m ^ one should expand the adiabatic 
potential En{p) (17) around po given by (23) using the extremal Ugxt up to the terms 
~ (p — Po)^- In such a way each adiabatic potential En{p) gives rise to the whole family 
of excitations in the quark subsystem. 

Strictly speaking, the adiabatic potentials (17) are valid either for n 1, or for 
p ^ where the spectrum (20) is exact. It is known, however, that quasiclassical 
results are usually rather accurate for small excitation numbers. Using with a lot of 
cautions the equations (22), (23) for the ground state energy Eo{p), we arrive to the 
effective one-dimensional adiabatic potential: 

= 21^5 +^;oW, (24) 

where = 15/4 defines the " hypercentrifugal" barrier for the lowest hyperspherical 
harmonic [10]. In accordance with equations (17), (18) adiabatic potentials diverge at 
small p , or at large a^, but this divergence is only the logarithmic one, as it is easily 
seen from the asymptotical behaviour of J(a) (see equation (19)). So for "realistic" values 
of quark masses hypercentrifugal barrier dominates the effective potential (24) at small 
values of p. Nevertheless, the spectrum of the ground state family differs substantially 
from the one given by the potential model (16), and it happens due to the account of 
junction zero oscillations. 

The confining force is not the whole story, and should be supplied by the Couloumb 
pair-wise interaction at small distances. Moreover, at small p the area law (7) is strongly 
contaminated by the effects of finite correlation length Tg, so that the string regime does 
not develop itself at full scale. The explicit calculations which take into account both 
these effects are in progress now, and will be reported elsewhere. It would be rather 
interesting to study the qqq system of light quarks in the presented approach, but this 
problem is much more formidable because the adiabatic approximation is not applicable 
for the light quarks. 

We acknowledge the extremely useful discussions with K.G.Boreskov, O.V.Kancheli 
and Yu.A.Simonov. This research is supported by the Russian Fundamental Research 
Foundation, Grant No 93-02-14937, and by INTAS - 93-0079. 
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